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Abstract 

The Buchdahl limit for static spherically symmetric isotropic stars is 
generalised to the case of five dimensional Gauss-Bonnet gravity. Our re¬ 
sult depends on the sign of the Gauss-Bonnet coupling constant a. When 
a > 0, we find, unlike in general relativity, that the bound is dependent 
on the stellar structure, in particular the central energy density and we 
find that stable stellar structures can exist arbitrarily close to the black 
hole horizon. Thus stable stars can exist with extra mass in this theory 
compared to five dimensional general relativity. For a < 0 it is found that 
the Buchdahl bound is more restrictive than the general relativistic case. 


1 Introduction 


An important question in general relativity is determining bounds on the mass 
and radius of stellar structures. The famous Buchdahl theorem [1] says that 
if we have a static perfect fluid solution to Einstein’s equation, whose energy 
density in non-increasing outwards, then the bound 


2M 8 



( 1 . 1 ) 


holds, where M is the mass of the fluid and R is its radius in Schwarzschild 
coordinates; defined by the location of the vanishing pressure surface. This has 
a number of important implications. For example, it tells us that the surface 
redshift is bounded, and that the boundary of a star always occurs after the 
apparent horizon in the Schwarzschild metric at r = 2M. 

Buchdahl’s theorem has been extended in various ways. It has been gen¬ 
eralised to include both charge and a cosmological constant, see for exam¬ 
ple Hill] and references within. Bounds have also been considered without 
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assuming Buchdahl’s assumption on non-increasing energy density [Slin]- Such 
inequalities have also been studied in modified theories of gravity, recently in [7] 
it was generalised to the case of f{R) gravity. Bounds have also been consid¬ 
ered in the context of Brane world scenarios HIS]. Studying mass radius bounds 
in modified theories of gravity is a way of testing their validity. If the theory 
predicts a bound which is too small or large one can experimentally measure 
whether there are stars violating these bounds. 

Bounds on the mass radius ratio in higher dimensional general relativity 
have also been considered by various authors. In d-dimensions the natural ratio 
to consider is 


2M 

Rd-3 


( 1 . 2 ) 


since this is the component appearing in the higher dimensional Schwarzschild 
metric, and thus allows one to derive bounds on quantities such as the stars grav¬ 
itational redshift [TU]. In [TT] Buchdahl’s theorem was extended to d-dimensions, 
where d > 4 and this was generalised to include a non-zero cosmological constant 
in |12] . Bounds were also considered without assuming Buchdahl’s asuumptions 
in d-dimensional spacetimes in m- 

Gravity in higher dimensions can be extended further than simply consider¬ 
ing higher dimensional general relativity. Gauss-Bonnet gravity is a particular 
natural theory to consider, and appears in the low energy effective action of 
string theory. This theory is a generalisation of Einstein gravity that adds an 
extra term to the standard Einstein-Hilbert action, which is quadratic in the 
Riemann tensor. When varying this extra term with respect to the metric only 
second order derivatives remain in the field equations, with the higher derivative 
terms cancelling out exactly, and thus the theory shares many of the nice prop¬ 
erties of general relativity. In four dimensions Gauss-Bonnet gravity and general 
relativity are equivalent, since the Gauss-Bonnet term in the action reduces to a 
total dervative, giving a surface integral and thus does not add a contribution to 
Einstein’s equation. But when analysing gravity in higher dimensions this extra 
term is non-trivial and it is thus natural to consider this extra Gauss-Bonnet 
contribution when considering higher dimensional theories. 

Many authors have considered fluid solutions in the context of Gauss-Bonnet 
gravity. Gonstant density interior solutions were investigated in |131114] . In [15] 
Boson star solutions were considered. Spherical symmetric gravitational collapse 
has been considered by many authors, see for example[Tni [U] and references 
therein. 

It has recently been claimed that Buchdahl’s theorem is not valid in Gauss- 
Bonnet gravity [18] . In this paper we investigate Buchdahl’s inequality in five 
dimensional Gauss-Bonnet gravity. Without the Gauss-Bonnet term the five 
dimensional Buchdahl inequality derived in m reads 


2M 3 
- 4' 


(1.3) 


It is of interest to investigate whether this Buchdahl bound, which uniformly 
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bounds the mass radius ratio away from the black hole bound, is a generic re¬ 
sult of gravitational theories, or is it in some sense a special property of general 
relativity. In this paper we use Buchdahl’s method [T] to show that a version of 
Buchdahl’s inequality does hold in Gauss-Bonnet gravity. In general relativity 
the interior solution saturating the Buchdahl bound is given by the constant 
density solution. This is not always the case for the inequality we have derived 
for the Gauss-Bonnet case. However it is found that for a stable stellar configu¬ 
ration, as in general relativity, the apparent horizon lies strictly below the radius 
of the star, however one cannot bound this uniformly away from the horizon. 

This paper is organised as follows. In Section II we review Gauss-Bonnet 
gravity and the five dimensional black hole solution and we explore the conse¬ 
quences on the mass radius bound for the constant density solution. In section 
III we derive the field equations for a spherically symmetric static perfect fluid 
and rewrite them in Buchdahl variables. In section IV we make the standard 
Buchdahl assumptions and derive bounds on the exterior vacuum metric. Fi¬ 
nally in section V we use these bounds to derive the corresponding Gauss-Bonnet 
mass radius bound. 


2 Gauss-Bonnet gravity 

In this section we introduce the action and field equations of Gauss-Bonnet 
gravity. The action for five dimensional Gauss-Bonnet gravity is given as follows 

S' = y [R + aLcB]^ + Sniatter- (^■^) 

The first term is the usual Einstein Hilbert action, whereas the last is the usual 
matter Lagrangian. The Gauss-Bonnet Lagrangian is given by the following 
particular combination of Ricci scalar, Ricci tensor and Riemann tensor 

Lgb = R^- AR^^Rab + RabcdR^^'^'^. (2.2) 

This particular Lagrangian appears from the low energy limit of heterotic su¬ 
perstring theory [19] . Note that we will not consider the effect of a cosmological 
constant in this paper. 

Varying the action with respect to the metric gives the following generalisa¬ 
tion of Einstein’s field equations 

Gab + oiHab = kTab- (2.3) 


Here Gab is the Einstein tensor, Tab is the usual energy momentum tensor, 
and the tensor Hab is given by 

Hab = 2 [RRab — 2RacRb ~ ‘2R^^Rabcd + Ra^^Rbcde] — -zgABLcB- 

(2.4) 
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From the interpretation of the low energy effective action from string theory, 
the coupling constant a is related to the inverse string tension which is positive 
definite, and thus the condition a > 0 is usually considered. This is also required 
for stability of Minkowski space in this theory. Nonetheless, in this paper we 
will also consider the case of a < 0, which has been considered by some authors, 
see for example HU HI]. We will also assume geometric units n = \ from now 
on. 


2.1 Gauss-Bonnet Black hole 


Assuming a static spherically symmetric metric and solving the field equations 
gives the following five dimensional black hole solution, first derived in [20] 


dr^ 


ds^ = -F(r)dt^ H-— + + sin^ 6 (dip^ + sin^ ip dijj^)) , 

F{r) 

where the function F(r) is given by 


(2.5) 


^ /) 16aM \ 




( 2 . 6 ) 


where M is related to the mass of the black hole or interior. Taking the limit 
a —^ 0 in this expression we recover the five dimensional Schwarzschild metric 


F(r) = l-^. 


(2.7) 


For this solution to describe the exterior of a star we require the apparent 
horizon to exist before the boundary of the star at r = R. The event horizon of 
a Gauss-Bonnet black hole is located at 


R = V2M - 2a 


( 2 . 8 ) 


This gives the following condition on the mass radius ratio 


2M 2a 


(2.9) 


assuming a > —R^jA. If a < —i?^/4 then we simply have the condition 


2M 1 

'W-2 


( 2 . 10 ) 


For positive a this mass-radius bound is less strict than the pure general rela¬ 
tivistic case. 

Birkhoff’s theorem in five dimensional Gauss-Bonnet gravity does not hold 
due to the presence of branch cuts, and thus this exterior solution is not unique. 
However we will only consider interior solutions which match to this exterior 
solution in this paper, as this is the branch which is asymptotically flat and 
agrees with the Schwarzschild metric in the limit a —>■ 0. 
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2.2 Constant density solution 

In standard general relativity the constant density solution saturates the Buch- 
dahl bound in any dimension. The constant density interior solution for five 
dimensional Gauss-Bonnet gravity was derived in |13| . We review this solution 
here and discuss the implications on the mass radius ratio. 

The energy density of the solution was taken to be a constant p and the 
pressure in the interior is given by 




1 - 


A* 


1-b 


2AVa 


B-y/ r2(l-/j)+4a 


where the quantity p is defined by 


p = \/l -I- Wawb, Wb 



( 2 . 11 ) 


( 2 . 12 ) 


The constants A and B and the mass of the fluid are determined by the matching 
conditions; requiring that the metric components gtt and prr are continuous at 
r = i?, as well ass gu,r being continuous on this surface. This implies that the 
mass is given by 


M = -pR^ 


and A and B are given by 


A = {l-B)^F{K), 


(2.13) 


(2.14) 

(2.15) 


In general relativity finiteness of central pressure of the constant density 
solution gives the Buchdahl bound (II.3|) . Let us examine what this condition 
gives in the Gauss-Bonnet case. Inserting r = 0 into (12.111) and requiring this 
to be finite and positive gives the inequality 

A 

- < -1. (2.16) 


This allows us to derive the following condition on F 

FiR) > ^ 

4(1 -I- ‘iawbY 


(2.17) 


The right hand side of this is well defined even in the case of a < 0, since 
1 -|-4auif, > 0. We will discuss the implications of this bound on the mass radius 
ratio in Section [5l 
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3 Field Equations for static star 


In this section we will consider the solutions for a more general perfect fluid. 
We assume a spherically symmetric static metric of the form 

ds^ =—e''di^ + e^dr"^ + (3.1) 

where is the metric of a 3-sphere and the functions n and A depend only 
on the radial coordinate r. We take the energy momentum tensor to be that of 
a perfect fluid 

= dia.g{-p,p,p,p,p) (3.2) 


where p is the energy density of the fluid and p is the isotropic pressure. The 
additional component of the energy momentum tensor in five dimensions, 
is required to be identical to and Tg under the assumption of spherical 
symmetry. Working in natural units with k = 1, the (t,t) component of the 
field equations ()2.3p gives 


P = ^ K - 2(1 - e^)) - ^(1 - e^) 


while the (r, r) component reads 
3e-^ 


p = 


2r2 


/ I 12ai/'e 

{riy' + 2(1 - e^))-3-(1 - e^). 


(3.3) 


(3.4) 


Conservation of the energy momentum tensor will give the following equation 


p' = -{P+P)v' 


(3.5) 


which is unmodified from general relativity. If these three equations are satisfied, 
the remaining Einstein equations, that is the 9 — 9, — and <p — p components, 

are identically satisfied, as in the standard four dimensional case. 

Let us write e~^ = 1 — f{r). From equation (13.31) we get 

{2a f + r'^f)' = ‘^p{r)r'^ (3.6) 

and hence integrating this gives 

2a f +r‘^f =‘^ J p{r'y^dr'. (3.7) 

Let us introduce the mass function 

p{r')r'^dr', (3.8) 

where the factor of 1/3 at the front required as the higher dimension mass 
function possesses an additional factor of l/{d — 2) [TT]. By matching this 
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interior solution with the exterior vacuum metric (1^ we see that evaluating 
this mass function at the boundary of the star gives the total mass of the fluid 

M = m{R). (3.9) 

Now we can solve for / to find the metric function e~^ 


e 


-A 


1 - / = 1 - 



16 am{r) 


4a 



(3.10) 


Now let us define our Buchdahl variables. First we introduce the function 
w{r) such that 

e~^ = 1 — 2r'^w{r) (3-11) 


which we can solve for w to give 


w{r) = 


\f^ 


16o'm{r) 


- 1 


8a 


We also introduce the following further variables x, y and 


C = e' 


V2 


y^ = e-^ = l- f{r). 


(3.12) 


(3.13) 


In terms of these new variables we can then rewrite the (r, r) component of 
Einstein’s equation (1331) as follows 


P = 6?/^^ — 6w + 48ay^w^. 


(3.14) 


Rewriting equation (13.61) in terms of w instead of / yields the following equation 
p = 6{xw^x + 2ic + 8aw{x)(xw^x + w)). (3.15) 

Inserting this into the conservation equation (1531) gives 


P,x = —{^{xw^x + ‘2-w + 8aw{xw^x + w)) + p)- 


c 


(3.16) 


By differentiating p.l4p with respect to x and inserting this into (I3.16P we 
can eliminate the pressure from these equations, which after simplihcation gives 
the following equation 


((1 + 8aw)yC,x),x - = 0. 

y 

Let us now introduce the final Buchdahl variable ^ given implicitly by 

dx 


d^ = 


y 


We can now use this to rewrite (I3.17|) as the following 

((1 + 8au;)C,e),{ - w^x^ = 0. 


(3.17) 

(3.18) 

(3.19) 
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4 Buchdahl bounds 


We are now in a position to derive our main results. In this section we will 
derive bounds on the metric function F at the boundary r = R. We will assume 
the standard Buchdahl assumption that the energy density is a monotonically 
decreasing function with respect to r. Thus we see from differentiating (13.121) 
that 


w'(r) 


rm'{r) — 4m(r) 

c; /i I 16 am(r) 

- 


< 0 


(4.1) 


In what follows we will denote quantities evaluated at the centre r = 0 with a 
subscript c and quantities at the boundary r — R with a subscript b. 

Now from inequality 63D we have 


tc.x < 0. (4.2) 

Inserting this into (13.191) yields the following relation 

((1 + 8au;)C,{),j = w.xC < 0. (4.3) 

Thus we can deduce that 


(1 + 8aw)C^^ > (1 + 8awb)iC,i)b- 


(4.4) 


Now evaluating at the boundary, using that ( is at the boundary so 
matches the derivative of the exterior solution, yields 


C« l&= 


F'iR) 


Wb 


AR 1 + 8awb 
Inserting this back into (|4.4I) thus gives 

(1 + 8aw)(^^ > Wb 

Integrating both sides of this with respect to ^ gives 

rib 


’ r 

(1 + 8aw)( f > 2wb / —dr 
Jo y 


(4.5) 


(4.6) 


(4.7) 


Now we find a lower bound for the left hand side of (I4J1) . Using the fact w is 
decreasing means y < \/l — 2wbr^ and hence 


■ dr > 


lo y 


r '' * 

Jo Vl - ^Wbr^ 

(4.8) 

7^(1 - \/l - 2wbr‘^) 

2wb 

(4.9) 


(4.10) 


Now to examine the left hand side of (14.71) we will need to consider the cases of 
positive and negative a separately. 
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4.1 Bounds fora > 0 


First we will consider the more physically relevant case a > 0. Now since ic is a 
decreasing function, (14.41) is a weaker inequality than the corresponding one in 
general relativity; and it will no longer necessarily be saturated by the constant 
density solution. 

Now we wish to find an upper bound for the left hand side of (14.7|) . Since a 
is positive and w is decreasing, w can be bounded by its central value 



rib 

{1 + 8aw)C,^^d^ < / (1 + 8aWc)C,cd^ 

Jo 

< (1 + 8awc)Cb 

= il + 8aw,)^/F{R) 


(4.11) 

(4.12) 

(4.13) 


And hence going back to (14.71) and using (14.131) and (14.101) we find the following 
inequality holds 


VF{R)il + 8awc) > (1 - VFiR))- (4-14) 


Rearranging this to find F(R) gives the following bound on the metric function 


F{R)> 


1 

4(1 + iawc)'^ 


(4.15) 


In the case of a constant density solution, we have that Wc = Wb, and hence 
our inequality (14.151) will in this case agree with the bound (j2.17l) found in the 
constant density case. 


4.2 Bounds for a < 0 


Now we consider the case a < 0. We must impose the condition 1 + 1^^ > 0 
in order for the exterior solution to be well defined, which implies 1 + 8awb > 0. 
Now, since w is decreasing, this time we can bound the left hand side of (I4A1) 
using the boundary value of w 

j {1 + 8aw) < {1 + 8awb) f C,id^ (4-16) 

Jo Jo 

< {l + 8awb)VF{R) (4.17) 

A lower bound on the right hand side of K7\i was already obtained, and thus 
we find 

(1 + 8awb)^yF{R) > (1 - ^/F{R)). (4.18) 

Hence we can derive the following bound on the metric function F{R) 

1 


F(i?)> 


4(1 + Aawb) 


2 ■ 


(4.19) 


This agrees with the bound (12.171) that was found for the constant density 
solution. 
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5 Mass-radius ratio bounds 


In this section we will examine what the implications of the bounds (14.151) , (I4.19|) 
on the metric function F(R) have on the mass radius ratio. 

First we will examine the case of a > 0, where the inequality (14.151) was 
derived. Let us define the quantity 

5a := 1 + AaWc > 1. (5.1) 


Now we can rewrite Wc in terms of the central energy density as follows 


Wr = lim 

r—vO 


1 + 


IQcxm 


- 1 


8a 


8a 


(5.2) 


Then rearranging (14.151) allows us to derive the following bound on the mass- 
radius ratio 


2M 1 2a 1 1 


(5.3) 


This is strictly weaker than the corresponding five dimensional general rela¬ 
tivistic Buchdahl bound. The right hand side of this inequality can in principle 
approach the horizon bound if the central energy density is large. However, this 
result is not necessarily saturated by any solution, as we have bounded the left 
hand side and the right hand side of (HTl) by different solutions. Thus in order 
to examine whether the bound does get arbitrarily close to the horizon bound, 
we will examine the constant density solution again. 

In the case of both the constant density solution for both positive and neg¬ 
ative a, and the generic case of a < 0, the same bounds (j4.19|) . ()2.17|) were 
derived 


F{R)> 


1 

4(1 -I- Aawb)^ 


(5.4) 


Now since the right hand side involves rcf,, which is given explicitly by 


Wb 



16Ma 

i?4 


-1 


8a 


(5.5) 


finding an inequality on the mass-radius ratio will reduce to solving a cubic 
equation. Rearranging (15.41) gives the following cubic equation in Wb which 
must be satisfied 


F[{wb) := —128a^R^wl + 64a (a — i?^) — 8 {R? — 4a) Wf, -I- 3 > 0 (5.6) 


Let us introduce the variables 

/ = ^R^Wb, 



(5.7) 
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In the limit a —>■ 0, / = 2M/R‘^. Now the inequality (15.61) becomes 
H{f) = + 16/3(/? - l)f + 4(4/3 - 1)/ + 3 > 0 


(5.8) 


For positive /3 the cubic has a positive discriminant and hence has three real 
roots. Only one of these roots is positive. One can see this by noting that the 
cubic is negative, H{—1/2(3) = — 1 < 0 and H{0) — 3 > 0 and hence two of the 
roots must be negative. Let us call this positive root 7 . Now since the horizon 
occurs and / = 1, and iJ(l) = —l< 0 we readily see that we must have 7 < 1 , 
and thus the bound on / lies strictly below the horizon. 

For negative (3, the discriminant of this cubic is negative when approximately 
— 1.32343 < /3 < 0. In this range the cubic has only one real root, which is 
positive, since the cubic is negative and iL(0) = 3 > 0. We will denote this root 
by 7 also. For (3 < —1.32343 the cubic has three real roots, all of which are 
positive. However, we note for negative (3 the inequality 


1 + 2/3/ > 1 + 4/3/ 



16aM 


> 0 


(5.9) 


holds, and thus / < —1/2/3. Now, H{—1/2(3) = —1 < 0, H'{—1/2(3) = 0, 
and thus the cubic has a turning point at this value. By looking at the second 
deivative of H we see that this turning point is a minimum for (3 < —1.32343 
meaning that two of the roots of the cubic are above —1/2(3. Thus only one 
of the roots of the cubic is less than —1/2(3, and we will also denote this root 
again by 7 . Hence for all values of (3 the cubic inequality reduces to the simple 
condition / < 7 . 

We plot the graph of this positive root 7 of this cubic in Fig[T] We note that 
7 varies continuously with /3; the root 7 belongs to the same branch cut. We see 
that for a > 0 the result rapidly approaches the horizon bound at Wb = 1. For 
negative a the bound on 7 is much stricter than the general relativistic case. 

Now this root of the cubic 7 does not have a pleasing analytic form. However, 
in terms of 7 we can find the following bound on the mass-radius ratio 


2M _ 2a7, 

-^ < 7( + -^) 


for B/^ > — 4 ( 37 , whereas for E? < — 4(37 have 


2M 7 

'W -2 


(5.10) 


(5.11) 


This is less than the horizon bound, but 7 —>■ 1 as a/B/^ gets large, and therefore 
there is no uniform bound below the horizon for (3 > 0 . 

But for a < 0 we can find some simple results on the mass radius bound. 
To do this we note a sufficient condition for (15.41) to be true is that 


f(R) > i 


(5.12) 
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Figure 1: Plot of the root 7 against fi := a/R^. The top dashed line indicates the 
horizon, whereas the lower dashed line indicates the general relativistic bound. 
The point shows the location of vanishing a, agreeing with the five dimensional 
general relativistic result. Positive /3 gives a less stricter bound than in general 
relativity, and approaches the horizon as /3 increases, whereas negative /3 gives 
a stricter bound and approaches 0 as ,5 decreases. 


This is equivalent to noting that in this case 7 < 3/4. Analysing this weaker 
bound alone allows us to draw important conclusions. Rearranging this, we find 
the two following bounds on the mass radius ratio, dependent on the radius of 
the star. For B? < —3a we have 


2M 3 

'W-% 

and for R? > —3a we have 

2M 3 9 

- 4 8i?2 “ 


(5.13) 


(5.14) 


6 Discussion 


In this paper we have investigated the equivalent of the Buchdahl bound in five 
dimensional Gauss-Bonnet gravity. For positive coupling constant a we have 
derived the inequality (15.31) . Unlike in five dimensional general relativity, the 
bound is not independent of the stellar structure since it depends on the central 
energy density through the value of pc- However, one important consequence 
of our bound is that it is strictly less than the position of the event horizon, in 
particular 


2M , 2a 


( 6 . 1 ) 


This is because the quantity da is finite and greater than 1. This means that 
Gauss-Bonnet black holes cannot have a perfect fluid interior either. 

In principle the right hand side of the inequality (15.31) can approach the 
horizon bound m if the central density becomes large. Thus unlike in general 
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relativity we cannot give a generic bound on quantities such as the surface 
redshift. In five dimensional general relativity, the surface redshift z is bounded 
above by z < 1. However our result shows that the redshift bound will be 
dependent on the specific stellar structure, agreeing with the conclusion in El- 
Also, unlike in general relativity, the constant density solution does not saturate 
the bound (|4.15p . However, even analysing the constant density solution, we 
see that this can be arbitrarily close to the black hole bound, and so we can 
see that this black hole bound is saturated by a fluid solution in Gauss-Bonnet 
gravity. 

In the limit a —>■ 0 of this inequality we recover five dimensional general 
relativistic Buchdahl bound 


2M 3 
- 4' 


( 6 . 2 ) 


The Gauss-Bonnet inequality is less restrictive than this general relativistic 
bound, and thus we are able to conclude that the appearance of the Gauss- 
Bonnet term allows stable configurations of stars with more mass in a given 
radius than their general relativistic counterparts. 

On the other hand, considering negative a gives a completely different 
situation. In this case these the constant density solution does saturate the 
bound (|4.19l) on the function F{R). The weaker inequality (15.41) is also true, 
and this allows us to derive the explicit bounds (15.131) and (15.141) , and are valid 
depending on the size of the radius of the star. Both of these bounds are 
stricter than the general relativistic bound, and are independent of the exact 
stellar structure. Again, the bound (I5.14p correctly reduces to the standard five 
dimensional Buchdahl bound in the limit a —>■ 0. This means we can fit less 
mass into a given radius and while maintaining a stable stellar structure. These 
inequalities will also imply a smaller upper bound for the surface redshift. The 
right hand side of inequality (15.1111 approaches 0 as the ratio a/B? decreases, 
and thus as a or i? decreases the mass-radius ratio gets very small. 

Of further interest would be to consider the maximum mass of a neutron 
star in Gauss-Bonnet gravity. In general relativity, it is known that Mmax < 
3.2Mq [22], with Mq the solar mass, where certain physical assumptions were 
made on the structure of the neutron star interior. Such an analysis has been 
done for higher dimensional general relativity for one particular EoS in |23j . 
where it was found that in higher dimension the neutron star maximum mass 
violated the Schwarzschild bound, and thus did not in fact describe a neutron 
star. Thus a further investigation of the structure of higher dimensional neutron 
stars in both general relativity and Gauss-Bonnet gravity should take place. 
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